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In the light of recent muon (g — 2) result by the E821 experiment at BNL, the 
importance of and interest in the total hadronic cross sections in e+e — collisions 
and il(e+e~ — » hadrons) have been heightened. We report on the integrated 
data compilation of available data of e+e - — > hadrons, transformed to a unique 
and meaningful style that implements the pure QED radiative corrections, from 
the PPDS and HEPDATA databases by the COMPETE Group. In particular, 
the continuum data sets are extracted where the direct comparison of the parton 
models and different variants of the pQCD corrections can be made and which will 
be useful for the future refinements of a B (Q 2 ) and &qed{Q 2 )- 



1 Introduction 

In spite of the tremendous experimental and phcnomenological efforts that 
have been devoted to the study of the hadronic production in e+e - collisions, 
it is simply the state of arts that there is no truly integrated view on the 
experimental situation for the main observable such as the total cross section 
for the reaction e+e - — > hadrons and the ratio 

er(e+e~ — > hadrons) 
er(e+e - — > u+u-)q E d 

To be specific, we note that some experiments published only the data of the 
cross sections, a, while others did only the R data. The situation is more 
complicated as for the QED radiative effects to the data, because various 
data sets have been corrected by different ways and to different degrees. For 
example, in constructing R, two different forms of theoretical calculations had 
been used, i.e., one with cxqed(0) and another taking into account the accepted 
form of evolution of olqed{Q 2 ) then. Furthermore, some data were never 
published in numerical form and have been read by compilers from graphs. 
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Determination of the precision data for cr(e + e~ — > hadrons) and i? has 
been a continuing endeavor in particle physics as new precision experiments 
became feasible and has received renewed interests recently. Some of the re- 
cent experimental activities that revived the interest in the data of c(e + e~ — > 
hadrons) and R are: the precision electroweak data that allowed us to cal- 
culate the radiative corrections and extract the QED coupling constant; the 
anomalous magnetic moment of the muon (g — 2) which showed some 2.6a 
deviation from the standard model prediction initially compared to the new 
result of the E281 experiment at BNlEl; and the r decay data. In the analy- 
sis of all these experiments, the hadronic contribution is the major source of 
the ambiguity and errors. A newj-controversy of the the light-by- light scatter- 
ing contribution (as for the signp to the vacuum polarization, the dominant 
leading contribution to the hadronic contribution reinforces only the need of 
precise values of these data. 

The aim of this worlfil is to prepare a comprehensive compilation of a and 
R transformed, whenever possible, to a unique, and meaningful, style by imple- 
menting the pure QED radiative corrections. The second aim is to construct a 
procedure to extract "continuum data sets" where the direct comparison of the 
parton model with different variants of the pQCD corrections can be made. In 
selecting data for the compilation we use the expert assessments from various 
recent reviews and works dealing-with precise estimation of running olqed{Q 2 ) 
and a s (Q 2 ) (see, e.g., Reffl, Reffl, and Reffl). 

2 Data Selection and Normalizations 

All published data were extracted in numerical form from the RPDS RD i 
database and from the Reaction Database of HEPDATA system u. Certain 
data were excluded based on the following exclusion criteria: 

• All preliminary data were excluded; 

• Data obtained under incomplete kinematical conditions, and not extrap- 
olated to the complete kinematic region, were excluded; 

• Data in the form of dense energy scans are omitted if the authors pub- 
lished in addition the data averaged over wider energy bins. 

The data sets that survived the above exclusion criteria were subdivided 
into the following six categories: 

<Tj d : Cross section data corrected in the original works for the contribu- 
tions of two-photon exchange diagrams, the initial state bremsstrahlung and 
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for the initial state vertex loops, but not corrected for the vacuum polarization 
contribution to the running chqed- These "semi-dressed" cross sections take 
into the full j /Z propagator with vacuum polarization effects. 

cr** 1 : Cross section data radiatively corrected according to the proce- 
dure of Bonneau and MartinEI in the original works. This procedure included 
radiative corrections for the initial state radiation, electronic vertex correc- 
tion, and the correction for the electronic loop in the photon propagator. 
Thus it partially took into account the vacuum polarization effects. To ob- 
tain from these cross sections the "semi-dressed" ones we rescale <t| by a 
factor 1/(1- Aa<Q ED (s)) 2 . 

R^ ave : Data on R obtained from measurements of the "semi-dressed" 
(o-2 d ) cross-section divided by the point-like muonic cross section with fixed 
a QED = &qed{Q) m the works. To obtain the "bare" i?-parameter we rescale 
these data by the factor (a QED (0) / 'a QED (s)) 2 / '(1 - Aa e QED {s)) 2 . 

R^ ave : Data on R obtained from measurements of the "semi-dressed" 
{c^ d ) cross section divided by the point-like muonic cross section, in the orig- 
inal works, but with the running o.qed(s) which thus obtained the "bare" 
i?-ratio. 

tr^ d : LEP I cross section data at the Z peak not corrected for initial 
state radiation and electronic QED vertex loops, and LEP II - III cross section 
data at -Js > 130 GeV with a cut s'/s = 1 - (Vs/2)(1 - 0.7225), Vs 7 being 
an effective mass of the propagator. These data were rescaled to the "semi- 
dressed" cross section cr^ d as follows. First, the theoretical cross sections^ 
cr*^ and u^om were calculated using ZFITTER 6.30 package E3, assuming 
the values of the standard model parameters from PDGllIl. Here <r'^ t means 
the cross section measured with cuts applied in the particular experiment and 
determined from the fit by ZFITTER and cr% L orn denotes the theoretical cross 
section, calculated by ZFITTER in the Improved Born Approximation (IBA) , 
that includes the full 7 / 1 Z propagator with vacuum polarization effects. Then 
the experimental cross section was multiplied by the factor cr^ l orn /o'*^ t . 

<Tg d : Low energy data {2m 7T < E cm < 2 GeV), the treatment of which 
requires special consideration. 

All the cross sections a(e + e~ — > hadrons) in this range are obtained by 
summing the exclusive channels and therefore give only a lower estimate of 
the total hadronic cross section. Below 1 GeV we summed up the 2ir and 3ir 

^Values for cr™. were calculated using ZFITTER subroutine ZUTHSM with argument set- 
tings according to the cuts applied in the LEP I-II-III measurements. Values for crfr 1 were 
calculated by the same subroutine but with the special flag switching it to calculate the IBA 
cross sections. 
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channels using a linear interpolation of the individual data sets within specific 
energy regions and combined them to give the total hadron cross sections in 
these regions. The errors were calculated according to these interpolation and 
summation procedures. 

Such an approach works well if all the data are evenly distributed over 
\fs and have comparable errors. Otherwise it may cause false resonance-like 
structures in the total cross section in the given y/s interval if there are few 
points of the leading channel with large errors, and if this interval is filled more 
densely by the points of the minor channels. Such false structures are possible 
and likely to occur for the exclusive channel data in the range 1.4 < y/s < 2 
GeV, where we summed the contributions from the channels > Shadrons, 
7r + 7r~ , K + K~ , and KgK^. To avoid such false structure we have determined 
the sum of exclusive channels not for all points, where at least one channel is 
measured, but in a fewer number of points, for which the channels yielding the 
major contributions are identified. In other aspects the summation procedure 
was the same as the one used for the data below 1 GeV. The exclusive sum data 
in the 0.81 < y/s < 1.4 GeV region are taken from the papeiS. The data used 
for the exclusive channel summation around 4> resonance (0.997 < y/s < 1.028 
GeV) were corrected for the initial state radiation, electronic vertex loops and 
leptonic 1-loop insertions into the 7 propagator in these data. We have prop- 
erly rescaled these data. In the remaining 0.81-^1.4 GeV data they applied 
QED-corrections using the Bonncau and Martin □ prescription and all these 
data points were rescaled as in the <r| d -case. 



3 Data Compilations on a sd and R bare 

After selecting and rescaling the data, normalized to the full 7/Z contributions, 
we assemble the complete data set of the total hadronic cross sections <r sd , in 
accordance with the symbolic relation 



/ rybare , , r>bare\ —t 1 ^ 

\.i u 4 ) ' ° QED pole, running a QED 



and the complete data set for the i?-ratios 

cr(e+e~ — > hadrons) 



j^bare 



er(e+e — > fi+fi )qed pole, 



running olqed 



As described above, in creating compilations which are uniform in the 
sense of implementing a standardized pure QED radiative corrections to the 



"raw" published experimental data, the values of the running aQED(s) and 
ActQ ED (s) are needed. The numerical values of these parameters and estimates 
for their uncertainties are described in the next section. Both compilations are 
stored and maintained in the PPDS CS database, corresponding to the record 
of the original data but with the rescaled data points. Also stored are brief 
descriptions of the applied conversion procedures in the special comment in 
each record. The er and R bare compilations obtained by this procedure are 
shown in the Figs. 1 and 2, respectively. 

4 chqedi ^ a QED an< ^ the Continuum Region 

4-1 Theoretical Relations. 

The QED running coupling constant oiqed{s) can be expressed in the form 



q(Q) 

aQED{s) = 1 - Aa( S ) ' 

Aa(i) = Aa had (s) + Aa lep (s) 



(1) 



where Aa had (s), Aa lep (s) are hadronic and leptonic contributions to the QED 
vacuum polarization, respectively. 

Aa lep (s) is well approximated in perturbative QED by the sum of leptonic 
loop contributions, 



Aa lep (s) = J2 



a(0) 

~3~r 



m — o — - 

to? 3 



s 



(2) 



The situation with respect to Aa had (s) is more complicated due to an essen- 
tially non-perturbative nature of the strong interaction at low energies. Using 
the unitarity and the analyticity of the scattering amplitudes, one can express 
Aa had (s) in the form of a subtracted dispersion relation 



Aa' lad (s) = - 



a(0)s f°° R(s')ds' 

3?r J 4m 2 s '( s ' - S-iO) 



(3) 



where R(s) is the "bare" hadronic ii-ratio. 

This relation allows to effectively evaluate the pQCD i?-ratios in the con- 
tinuum y/s intervals in terms of experimental R data containing the non- 
perturbative effects. 
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We have evaluated the dispersion integral in two parts, 



Aa had (s) 



a(0)s 
3tt 



ds ' 4(4 

Ami. s \° 



Tidat: 
^bare 



lis') 



iO) 



ds 

Scut 



s'(s' 



iO) 



(4) 

where the first integral was calculated numerically as a trapezoidal sum over 
the weighted average of experimental R points in the range 2m 7r < ^fP < 
^Tt = 19.5 GeV c 

Our numerical evaluation procedure is in general similar to the one applied 
inRef.0.0 

We obtained Aa /lQd (Af|) = 0.0274±0.0004(exp.), which is consistent with 
the results 0.027382 ± 0.000197 and 0.027612 ± 01)00220, obtained with two 
methodically different low energy data sets in Rcf.ta. This result is preliminary 
and serves only as a guarantee against egregious errors in our procedure. The 
procedure of Aa had (M z ) evaluation requires further refinement. 

In the parton model, before QCD corrections are applied, the i?-ratio is 
given by 

1 



R 



3> R°„ 



Q 



P q {l + ^{l-0 l q ))-R v q v +PtRf A 



(5) 



where 



with 



R 
R 



vv 

<? 

AA 
1 



2 2 
it' 



+ 2e e e q v e v q Rex + (v e + a e )v \x\ 



(G) 



X 



16s 2 c 2 {s-M 2 z + iM z T z ) 



c Choosing s cu t we must satisfy at least the following two requirements: 

1) s cu t should lie well above all the hadronic resonances, i.e. 1-loop perturbativc QED can 
be used to calculate hadronic vacuum polarization at s > s cu t\ 

2) there should be enough densely distributed experimental points at s < s cu t for the 
trapezoidal numerical evaluation of the dispersion integral over 2m n < < y?s cu t- As 
there are few experimental points in the interval 13 < \fs < 30 GeV, yJs C ut = 19.5 GeV 
appears to be a compromise between these requirements. 

d As the calculation of fi|J°'° from the original a and R data (see Section 1) in turn requires 
the evaluation of Aa had (s), we applied the following simple iterative procedure: Taking as 
the zeroth approximation R(s) as obtained from the original data using 1-loop perturbativc 
Q.qed{s), we evaluate Aa h ° l! ' 1 '(s), which in turn is used to compute the next approximation 
of R(s) from the original data, and then evaluate the integral again to obtain Aa had ^(s), 
and so on. This process converges well and after the two iterations the algorithmic error of 
Aa (s) in the continuum region is much less than the error coming from the experimental 
data uncertainty. 
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v = VP (T 3 - 2Q1 2 ) 
a= ^pT 3 

P = 1 + 1R 2 TO t + • • • ■ (7) 

Here s 2 , c 2 are effective sin 2 9w , cos 2 6*vf defined through renormalized cou- 
plings at s — M§ and 

V2G_f ■ p ■ s ■ c 

The dominant correction to p originates from f-quark loops in W and Z 
propagator due to the large mass splitting between the b and t quarks which 
represents an SU(2) breaking in the unified electroweak. theory. 

Including QCD loops, R is given by (see, e.g., RcfE£l and RefEZI): 



R = 3 £ R, 
= 3£ 



R o + (^. + ^ (1.9857-0.1153^/) 
(— ) 3 (6.6368 + 1.20017V/ + 0.00527V 2 + 1.2395(E Qqf) 



{hR v q v + f 2 R* A ) 



(9) 



where Nj is the number of active quark flavors. The coefficients fx and /2 
depend only on the quark velocity (3 = (1 — 4M 2 /s) 1 / 2 —where M q is the 
threshold for the flavor production of quark q. Schwingerll3 calculated fi in 
the QED context, which is accurate enough for our purposes: 



h 



The coefficient / 2 (which have no counterpart in QED) were calculated by 
Jersak et alE3, which we parameterized as 

f 2 =a A p A q +a 3 fi + (l-a 4 -a 3 )f3 2 q (11) 

with 04 = —16, a 3 = 17. From equations (||), (||), ^) one can easily see that 
the relative error of R is less than 10 -7 at ^/s = 10 GeV and less than 0.5% 
at the Z pole, even if we allow a 100% error for f 2 . As we did not use Z pole 
data in our fits, this could not be a source of large theoretical errors. 

In the massless quark limit, f± = f 2 = 1. However, at y/s = 35 GeV for 
b quark one has (3 = 0.963 and thus /i ~ 1.3 and f 2 — 1-7. Indeed, such a 
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parameterization of mass effects is valid in QCD only in 0(a s ) order | _ 1 Now 
the .correct parameterization is known up to O(a^) order (see, e.g., Reft3 and 
ReftD), but it has not been implemented yet into our programme. It may 
be that it results in large enough discrepancies just in the region where a s /ir 
becomes large. 

We used the following three- loop parameterization for a s (see, e.g., Rcfl 14 



In (Q 2 /A 



with 



MS' 



An 

Poa s 



I Pi . 



An 

Po^s 



tan 



Pi 
Pi 
Pi 
PI 



An 

PoOls 

2p2 

PI 



+ Pl 

Pi 

Poa s 

An 



(12) 



Po 
Pi 

P2 



11 - 2JV//3 , 
2(51 - 19JV//3) , 

(2857 - 50337V//9 + 325iV) /27) /2 , 



^Afa/Pl-Pllft 



The onset of a new flavor q takes place at Q 2 



(13) 

Am 2 , where, in general, 

M q . There are five different Ajj^(Nf) corresponding to Nf = 2, 3, 4, 5, 6 
above the appropriate quark thresholds. To sew up a s at the thresholds we 
need to use a matching condition for A's: a s (Q 2 ) is determined by fixing, say, 
a s (Mz), which is the parameter determined from the fits in our work. 

We solved this equation at given Q 2 numerically by rewriting it in the form 
x = f(x), where x = 47r/(/?oa s ). It enabled us to obtain correct a s even at 
Q 2 /A 2 - 1, where ihe expansion by powers of leading logarithms is no more 
valid (see, e.g., Rem). At Q 2 /A 2 ^ 1 both methods yield the same result (see 
Fig. 4). 



4-2 Preliminary Fit Results 

In order to check the consistency of the description of the data on 

R = a(e + e~ — * hadrons) I o(e + e~ 

with the theoretical _R-ratio expression we have performed several fits using a 
sub-set of the compiled data set of R which is applicable to the perturbative 
domain. This sub-set of the data is defined by the the following steps: 
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1. The low energy region 2m ff < s/s < 2 GeV is completely excluded; 

2. All the data above 70 GeV are excluded, as we did not attempt to perform 
any fits of Z pole parameters Mz, ^z, ~S\v an d M top - This exclusion is 
justified because the data above 70 GeV weakly depend on a s (Mz) and 
the quark masses and thus their influence on the fit of Z pole parameters 
would be rather negligible; 

3. Data covering the narrow hadronic resonances, J/ip(lS), ip(2S), 
V<3770), T(15), T(2S), T(3S), T(4S), T(10860) and T(11020) are com- 
pletely excluded. (The visible widths of these resonances are determined 
by the machine energy spread AEb eam ~ 10 MeV.); 

4. Also completely excluded are the SLAC-SPEAR-MARK-1 data in the 
2.6 GeV < y/s < 7.8 GeV range. These data systematically lie ~15% 
above other experiments in the same interval. 

We distinguished the 7* — > qq production threshold quark masses M q from 
the QCD quark "masses" m q determining the energy scale of the onset of the 
new flavour q in the running of as{s). After performing the exclusions 1) - 
4) we fixed a s (M z ), M z , T z , Syy, M q (q = u,d,s,t), m q (q = u,d,s,c,b,t) at 
the values shown in the Table 1, leaving just M c and Mf, as free parameters. 
Another essential feature of our fit (I) was that we retained as free parameters 
the left and right boundaries of the interval [ v /sl, -^/si], 3.0 < -/si < 3.670 
GeV, 3.870 < yfs2 < 5.0 GeV, which was to be excluded from the data set in 
the process of x 2 /dof minimization itself^]. Thus, the number of the degrees 
of freedom was also variable during this fit, for which the standard MINUIfO 
program is used. Such a x 2 /dof minimization procedure cancelled possible ar- 
bitrariness in the exclusions of broad resonances in the region above cc thresh- 
old. 

In the fit (II) we released a s (Mz) free, fixing at the same time the excluded 
region y/s — 3.09-^4.44 GeV. M c , Mb were fixed at their fitted values obtained 
in the fit (I). Fit (II) gives somewhat high value a s {M z ) = 0.128 ± 0.032. 



e Narrow ip family resonances between 3 and 4 GeV were excluded by default. Somewhat 
controversial is the question, whether the gaps between tp's can be really treated as the 
continuum regions. R ratio demonstrates no step up until yfs = 3.9 — 4.0 GeV, just at the 
left knee of the broad -0(4040) resonance, where the continuum approximation of R does not 
work. Nevertheless, this step results in the fitted c-quark mass M c = 1.971 (see fit (III) in 
the Table 1) 
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5 Summary. Assessed Data Compilations. 

In summary: 

• We have created two complementary computerized "raw" data compi- 
lations on a and R with data presented as in the original publications 
in all cases except the low energy subsample. In the low energy region, 
where there are no direct measurements of the total cross section, we 
obtain estimates of the total cross section either as the sum of exclusive 
channels or as the sum of the two-body exclusive channels with the data 
on e + e~ — >> 3 hadrons. 

• On the basis of above two data sets we have created compilations of data 
on cr 3d and R bare with one-to-one correspondence between the data 
points. Where necessary the data have been rescaled to the standard 
style of implementing the pure QED radiative corrections to the initial 
state and to the photon propagator to produce a data set which is suitable 
to test the parton model with pQCD corrections to R bare . and to be able 
to obtain more reliable estimates for AaQ^ D (Q 2 ). 

• From the total data compilation on R bare we have defined a "continuum" 
data compilation sub-set which can be used in conjunction with other 
data in the refinements of the a s {Q 2 ) evolution and possibly in the global 
fits of the Standard Model. 

All data files are accessible by: 

http://wwwppds.ihep.su:8001/comphp.htm] 
and will be accessible from the PDG site soon. 
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Table 1. Fit results in the defined "continuum" regions. Adjustable param- 
eters are shown in bold. 



Parameter 


I 


II 


III 




0.1181 


0.128(32) 


0.126(37) 


M z 


91.187 


91.187 


91.187 


T z 


2.4944 


2.4944 


2.4944 


s 2 


0.23117 


0.23117 


0.23117 


M t 


174.3 


174.3 


174.3 


M u 


0.140 


0.140 


0.140 


M d 


0.140 


0.140 


0.140 


M s 


0.492 


0.492 


0.492 


M c 


1.500(18) 


1.5 


1.9710(1) 


M b 


5.23 ± 1.16 


5.23 


6.0 ± 1.4 


m u 


0.003 


0.003 


0.003 


m d 


0.006 


0.006 


0.006 


m s 


n i on 


n 1 on 

U. 1ZU 


n i on 


m c 


1.2 


1.2 


1.2 


m b 


4.2 


4.2 


4.2 


m t 


174.3 


174.3 


174.3 




2m n + 2.0 


2m^ 2.0 


2m n + 2.0 




3.093-^3.113 




3.093^-3.113 
3.684-^3.688 
3.670-3.870 


Excluded 


3.175(247) -e- 
4.319(105) 


3.09 -=-4.44 


4.000^-4.400 




9.450-^9.470 


9.450^9.470 


9.450^-9.470 




10.000^-10.025 


10.000^-10.025 


10.000-e-10.025 


intervals 


10.34 -e-10.37 


10.34 -e-10.37 


10.34 -e-10.37 




10.52 -e-10.64 


10.52 -e-10.64 


10.52 -e-10.64 




10.75 -e- 10.97 


10.75 -=-10.97 


10.75 -e- 10.97 




11.00 -e-11.20 


11.00 -e-11.20 


11.00 -e-11.20 




70 -e- 188.7 


70 -e- 188.7 


70 -e- 188.7 


X 2 /dof 


0.690 


0.665 


0.822 
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Figure 1: World data on the total cross section of e+e - — > hadrons. The experimental data 
on J /if), ip(2S), and T(nS),n = 1..4 resonances is not shown. Curves are an educational 
guide. Solid curve is the prediction of the cross section in the three-loop QCD approximation 
with non-zero quark masses. Dashed curve is a "naive" quark parton model prediction. 
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Figure 2: i?-ratio. Data set is the same as in Figure 1. Solid curve is the ij-ratio prediction in 
the three-loop QCD approximation with non-zero quark masses. Dashed curve is a "naive" 
quark parton model prediction for the ratio parameter R. 
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Figure 3: Data set used for the preliminary fits. Curves are the same as in Figure 2. 
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Figure 4: Comparison of i? parameters obtained with a s (s), evaluated by our numerical 
method (solid curve) and by the method described in [5] (dashed curve). 
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